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ABSTRACT 

We obtain Keck HIRES spectroscopy of HVS5, one of the fastest unbound stars in the Milky Way 
halo. We show that HVS5 is a 3.62 ±0.11 M Q main sequence B star at a distance of 50 ± 5 kpc. The 
difference between its age and its flight time from the Galactic center is 105 ± 18(stat)±30(sys) Myr; 
flight times from locations elsewhere in the Galactic disk are similar. This 10 s yr 'arrival time' between 
formation and ejection is difficult to reconcile with any ejection scenario involving massive stars that 
live for only 10 yr. For comparison, we derive arrival times of 10 7 yr for two unbound runaway B 
stars, consistent with their disk origin where ejection results from a supernova in a binary system or 
dynamical interactions between massive stars in a dense star cluster. For HVS5, ejection during the 
first 10 7 yr of its lifetime is ruled out at the 3-tr level. Together with the 10 8 yr arrival times inferred 
for three other well-studied hypervelocity stars, these results are consistent with a Galactic center 
origin for the HVSs. If the HVSs were indeed ejected by the central black hole, then the Galactic 
center was forming stars ~200 Myr ago, and the progenitors of the HVSs took ~100 Myr to enter the 
black hole's loss cone. 

Subject headings: Galaxy: halo — Galaxy: center — Galaxy: kinematics and dynamics — stars: 
early-type — stars: individual (SDSS J091759.47+672238.35) 



1. INTRODUCTION 



Hills (|1988l ) first predicted unbound "hypervelocity" 
stars (HVSs) as the inevitable consequence of 3-body in- 
teractions close to the tidal radius of a massive black 
hole. There is overwhelming evidence for a 4 x 10 6 
central black hole in the Milky Way (|Ghez et all 120081 : 
iGillessen et al.ll2009l ). Theorists ex pect that the blac k 
hole ejects ~10" 4 HVSs yr^ 1 (e.g. iPerets et all l2007h ■ 
whic h means there are t housands of HVSs in the outer 
halo. iBrown et al.l (|2005l ) discovered the first HVS, a lu- 
minous B-type star traveling twice the Ga lactic escape 
veloci ty at a distance of ~100 kpc, and IBrown et al.l 
(|2012f) have subsequently discovered 15 more unbound 
B-type stars in their targeted HVS survey. Establish- 
ing the evolutionary state of the HVSs is important for 
establishing their ages, distances, and flight times. We 
define the difference between a HVS's age and its flight 
time as the 'arrival time' (t aTI ), the time between its for- 
mation and ejection. In this paper we derive i arr for both 
HVSs and unbound runaway stars. 

The arrival time provides a useful discriminant be- 
tween proposed ejection mechanisms. If HVSs are ejected 
in three-bo dy interact ions with the Milky Way's central 
black hole (Hills 1988), then the arrival times reflect the 
timescale for HVSs to achieve orbits that interact with 
the central black hole. For HVSs formed in the cen- 
tral region of the Galaxy, we expect t a rr= 0.1-1 Gyr 
(jMerritt fc Poonl 120041 : 1 Wang fc Merrittll2004| ). On the 
other hand, in both mechanisms for ejecting runaway 
stars from the Galactic disk - a supernova in a binary 
system or a dynamical interaction among massive stars 
in a dense star cluster - a maximum i arr ~ 10 Myr is set 
by the main sequence lifetime of >10 M Q stars. Thus, 
measuring i arr for an ensemble of HVSs should distin- 



guish between a Galactic center or Galactic disk origin. 

The evolutionar y state of most known HVSs 
i|Bio\vii el al.l [20 1 2 ) is ambiguous because their effective 
temperatures and surface gravities are consistent with 
both old, evolved stars (blue horizontal branch stars) 
and short-lived main sequence stars. Thus we must turn 
to other measures to establish their nature. Metallicity 
is one possibility; we expect that recently formed stars 
should have solar or super-solar metallicities. Metallicity 
is inconclusive, however, given the observed metallicity 
distribution function of stars in the Milky Way. 

Projected stellar rotation wsini is a better discrim- 
inant between evolved stars and main sequence stars. 
Blue horizontal branch stars have evolved through the 
giant branch phase and have median wsinz = 9 km s _1 ; 
the most e xtreme blu e horizontal branch star rotates at 
40 km s _1 ( Bchr 2003). Late B-type main sequence stars, 
on the other hand, have median v sin i = 150 k m s -1 ; the 
most extreme objects rotat e at >350 km s" 1 (lAbt et al" 



[20021: iHuang fc Giesl [20061 ) . lLockmann fcB aumgardt 



(2008) argue that HVSs may be spun-up by a binary 
black hole ejection, but there is presently no evidence 
for a binary black hole in the Galactic center. Close bi- 
naries that produce HVSs and runaways in the Milky 
Way may exhibi t slower stellar rotation because of tidal 
synchronization; Hansen] ()2007f l predicts that late B-type 
HVSs ejected by the Hills mechanism should have v sin i 
= 70 — 90 km s _1 . In any case, fast rotation is the sig- 
nature of a main sequence star. 

Of the B-type HVSs discovered to date, only 
HVS3, HVS7, and HVS8 have been studied with 
high-resolution spectroscopy. In all cases they 
are mai n sequence B stars with 55<fsinz<260 
km s" 1 (lEdelmann et all 120051: iPrzvbilla et all 12008b : 



km s" 1 (lEdelmann et~all 120051: iPrzvbilla et ail 12008b 
iBonanos et al.1 120081 : iLopez-Morales fc Bonanosl 12008 
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iPrzvbilla etaLl l2008d ). Moderate-d lspersion spec- 
troscopy of HVS1 s uggests it has vsini = 190 km s _1 
(jHeber et al.ll2008bl ). another short-lived B star. 

Here, we describe high resolution spectroscopy of 
HVS5, a g = 17.9 mag star located at declination +67° 
accessible only with Keck HIRES. HVS5 is a rapidly ro- 
tating 3.6 Mq main sequence B star. The difference be- 
tween its age and its flight time from the Milky Way is 
105 ± 18(stat)±30(sys) Myr, inconsistent with ejection 
models involving massive stars. 

In Section 2 we describe the observations and stellar 
atmosphere analysis. In Section 3 we discuss the arrival 
times for the HVSs and unbound disk runaways. We 
conclude in Section 4. 

2. DATA 

2.1. Observations 

We observed H VS5 using the HIRES spectrograph 
(|Vogt et all 119941 ) at the 10 m Keck 1 telescope, ob- 
taining four 1800 sec exposures the night of 2012 Jan 
29, and five more the following night. Both nights were 
clear with 0.7-0.8 arcsec seeing. The red HIRES colli- 
mator was used in an instrument configuration that gave 
spectral coverage from 3920 to 8350 A. A 1.1 arcsec slit 
gives a spectral resolution of R =34,000, with 6.7 pix- 
els per spectral resolution element. There are small gaps 
between the three CCDs that form the detector mosaic, 
sometimes resulting in the loss of all or part of a single 
echelle order. 

We also observed eight B stars selected from lAbt et al.l 
((2002T) that span a wide range in luminosity class (I, III, 
and V) and a large range in projected rotational veloc- 
ity (5 < wsini < 285 km s" 1 ). The stars are HR1328, 
HR1333, HR1399, HR1419, HR1420, HR1462, HR1573, 
HR1595, and HR1640. We used high S/N spectra of 
these eight B stars to provide a comparison for HVSS 
and to validate our analysis below. 

We used the pipeline package makeeQ to remove the 
instrumental signature, extract a one dimensional spec- 
trum for each echelle order, and calibrate the wavelength 
scale from Th-A arc spectra taken at the beginning and 
end of each night. Each exposure of HVS5 was individu- 
ally processed through makee, and the results summed. 
With a total integration of 4.5 hours, we achieved a S/N 
ratio of 70 per spectral resolution element at 4500 A. 

2.2. Spectral Analysis 

Visual inspection reveals that HVS5 is a fast rotator. 
The full width at half maximum of the Mg 11 A4481 
line compared with the lAbt et ahl (|2002| ) B stars sug- 
gests a projected rotation velocity of ~130 km s _1 . We 
turn to the most recent ATLAS9 model atmo sphere grids 
(|Castelli fc Kuruczll2004l ICastelli et al.lll997[ ) to perform 
a more quantitative analysis. W e use the program spec- 
trum (|Grav fc Corballvl I1994D to calculate synthetic 
spectra under the assumption of plane-parallel atmo- 
sphere and local thermodynamic equilibrium. We adopt 
macro- and microturbulunce velocities of and 2 km s , 
respectively. All synthetic spectra are smoothed to a res- 

1 MAKEE was developed by T.A. Barlow specifically for reduction 
of Keck HIRES data. It is freely available from the Keck HIRES 
home page www2.keck.hawaii.edu/inst/hires 
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Fig. 1. — Observed Mg II A4481 A line (upper panel) and its 
residuals (lower panel) compared to the best-fit tisini= 133 ±7 km 
s~ 1 model (red line). For reference, the green and blue lines are 
the ±3-<T models. 

olution of R = 34, 000 and re-sampled with spline inter- 
polation to match the sampling of the observed spectra. 

Our approach is to analyze the spectral lines on an 
order-by-order basis. We normalize the continuum, cal- 
culate the x 2 of each synthetic model against the data, 
and then fit the resulting distribution of x to derive 
the best-fitting parameters. Uncertainty estimates are 
provided by a sta ndard A% 2 offset to the minimum x 2 
([Press et al.1 119921 ) . Our final values are the weighted 
means and error-on-the-means of the parameters mea- 
sured from lines in different orders. 

We begin by using all of the spectral lines to solve for 
the heliocentric radial velocity. The best- fit +552 ± 3 km 
s _1 velocity is in perfect agreement with the +553 + 9 km 
s" 1 velocity measured from medium-reso lution spec- 
troscopy at the MMT (|Brown et al.l I2012H . The con- 
stancy in velocity is consistent with HVS5 being a single 
star, as one expects for the product of a binary disrup- 
tion. The radial velocity corresponds to a minimum ve- 
locity of +663 km s _1 in the Galactic rest frame (see 
iBrown et al1l20Tl . 

Next, we measure projected rotation using Mg II A4481 
A, the strongest metal line in the spectrum (see Figure[l}. 
The best-fit usini is 133 ± 7 km s . For comparison, 
Figure Q] plots the ±3-<r values as well as the residuals to 
the best- fit vsmi. The observed wsini is consistent with 
the m edian t>sin i of comparable B-type m ain sequence 
stars (|Abt et al.ll200l iHuang fe Giesll2006D . 

Given the observed wsini, we measure effective tem- 
perature and surface gravity by fitting the widths and 
depths of the T ff- and log (/-sensitive hydrogen Balmcr 
lines. We note that the model hydrogen lines are com- 
puted usiMtieD . Peterson routine adopted by SYN- 
the (|Kuruczlll993| ). which includes Stark and resonance 
broadening and fine structure in the cores. The best-fit 
values are T off = 12, 000 ± 350 K and log g= 3.89 ± 0.13; 
Figure [H compares the best- fit model with the data. 

Finally, we attempt to constrain the metallicity. Be- 
cause of the large usini, Fe lines are faint continuum 
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Fig. 2. — Observed hydrogen Balmcr lines compared to the best- 
fit model (red line) . The temperature- and surface gravity-sensitive 
lines give best-fit values of T oft = 12. 000 ±350 K and log 3= 3.89 ± 
0.13. 

fluctuations and are thus too weak to provide signifi- 
cant constraint. Si 11 lines at A4128, A4131, and A5056 
A are stronger and yield a best-fit Si abundance of 
[M/H]=-0.4 ± 0.5. The Mg n A448f A line (see Fig- 
ured]), on the other hand, yields a best-fit Mg abundance 
of [M/H]=+0.3 ± 0.5. Given the large uncertainties, we 
conclude that HVS5 is consistent with solar abundance. 

Figure [3] compares the measure d T e g and log g with the 
latest IGirardi et all (|2002t 120041) solar metallicity main 
sequence tracks. The ellipse in Figure [3] is the 68.3% (1- 
cr) confidence region. Interpolating the tracks indicates 
that HVS5 is a 3.62 ± 0.11 M Q star. As an illustration 
of the systematic un certainty, we derive 3.58 Mq from 
lEkstrom et al.l (120121) trac k s with rotation, and 3.72 M 
from IGirardi et all (|2002l . I2004D +0.2 dex super-solar 
tracks. These values are consistent within our 1-a un- 
certainty, thus the inferred mass is relatively insensitive 
to rotation and metallicity. 

Table [T] summarizes our stellar parameters for 
HVS5. We also list the parameters for HVS7 , 
HVS8, and HVS1 measure d bvlPrzvbilla et all (12008^ 
lLopez Morales fc Bonanosl (|2008l ). and iHeber et al.l 
2008 rj), respective ly. HVS3, a probable blue straggler 



( Bro wn et al.l [201 Oh. is not directly comparable and is 
not included in our discussion. 

3. DISCUSSION 
3.1. Hypervelocity Star Flight Times and Ages 

We use the IGirardi et al.l (|2002L 120041 ) tracks to derive 
age and luminosity from the measured T e ff and logg. We 
note that the T e s-logg error ellipse becomes a banana 
shape in the age-M 9 plot (Figure U]) because of the time 
evolution of these parameters. Interpolating the tracks 
indicates that HVS5 has an age of 170 ± 17 Myr and an 
absolute magnitude of M g = —0.65 ± 0.19. 

Knowing the luminosity of HVS5, we can calculate its 
distance and flight time. HVS5 has an apparent de- 
reddened magnitude of g = 17.557 ± 0.021 and thus a 
heliocentric distance of 44 ± 4 kpc. Assuming the Sun 
is located 8 kpc from the Galactic center, HVS5 has a 
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Fig. 3. — Measured T e ff , log g and the 1-a erro r ellipse for HVS5 
(in red) compared to Girardi ct al. (2002, 2004) solar metallicity 
main sequence tracks for 2.5-5 Mq stars (solid black lines); the 
tev = 75 Myr isochrone is plotted for reference. HVS5 is a 3.62 ± 
0.11 M star. 
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Fig. 4. — Same tracks as Figure [3] but now plotting age 
versus absolute magnitude. HVS5 has a formal age of 170 ± 
17(stat)±30(sys) Myr. 

Galactocentric distance of roc = 50 ± 5 kpc. We then 
take the Galactic potential model of lKenvon et al.l ([2008D 
and calculate all possible trajectories that HVS5 could 
have given its present distance and radial velocity. The 
trajectory that passes through the Galactic center has 
a flight time of ice = 65 ± 7 Myr. The error comes 
from propagating the distance and radial velocity errors 
through the trajectory calculation. 

Our Galactic center flight time estimate is appropriate 
for a wide range of Milky Way starting locations because 
HVS5 is located at high Galactic latitude and at large 
distance. Moving the assumed starting location of HVS5 
from roc = kpc to tqc = 10 kpc changes the flight 
times by ±8.5 Myr, which is similar to the estimated 
uncertainty in flight time. 
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TABLE 1 

Hypervelocity Star Properties 



HVS5 HVS7 HVS8 HVS1 



T eff (K) 12000 ± 350 12000 ± 500 11000 ± 1000 11000 ± 500 

logg (cgs) 3.89 ± 0.13 3.8 ± 0.1 3.75 ± 0.25 3.74 ± 0.20 

vsini (km s" 1 ) 133 ± 7 55 ± 2 260 ± 70 190 ± 40 

Mass (M e ) 3.62 ± 0.11 3.79 ± 0.09 3.49 ± 0.22 3.50 ± 0.18 

Age (Myr) 170 ± 17 170 ± 15 220 ± 25 220 ± 20 

M g (mag) -0.65 ± 0.19 -0.94 ± 0.14 -0.77 ± 0.35 -0.80 ± 0.29 

go (mag) 17.557 ± 0.021 17.637 ± 0.019 17.939 ± 0.016 19.688 ± 0.023 

roc (kpc) 50 ± 5 53 ± 4 60 ± 10 130 ± 18 

tec (Myr) 65 ± 7 105 ± 10 120 ± 20 175 ± 25 

i arr =Age-t GC (Myr) 105 ± 18 65 ± 18 100 ± 32 45 ± 32 



For self-consistency, we also derive the ages and flight 
times of HVS7, HVS8, and HVS1 using the same tracks 
and methodology. Table [T] summarizes the derived val- 
ues. The HVSs have ages of 170-220 Myr and flight times 
of 45-175 Myr. 

3.2. Links to Unbound Ejection Processes 

There are many ways to eject stars from their place 
of origin, but few processes can accelerate stars to un- 
boun d velocities . Be cause most B stars are binaries 
(e.g. IChini et"a^[20l2l) . disk "runaway" B stars are ex- 
plained by binary disruption mechanisms. In the case 
of a supernova in a binary system, the timescale of the 
process is the lifetime of the >10 M Q star that ex- 
plodes, 10 6 -10 7 yr. In the case of dynamical 3- and 
4-body encounters, e.g. in young star clusters, mas- 
sive stars are necessary to attain the unbound veloci- 
ties of HVSs and thus the timescale of the dynamical 
process is also ~10 7 yr. Except in rare circumstances 
(iGualandris fc Porteeies Zwartll2007t iGvaramadze et all 
120091 : ISilva fe Napiw otzki 2011), no runaway mechanism 
is expected to yield unbound velocities for 3 Mq stars 
(|Portegies Zwart|[2000t fPerets fe Subrll2012D . 

A more energetic and higher ejection rate process ex- 
ists in the Galactic cen ter: HVSs ejected by the central 
black hole (jHillsl 119 88) . The B stars that presently or- 
bit Sgr A* on short-period, eccentric orbits are, in this 
scenario, the former companions of HVSs; their progen- 
itors are believed to have formed further out and then 
moved in t owards the black hole through dynamical pro- 
cesse s (e.g. lPerets et al.ll2007l I2009T: iMadigan et al.|[200l 

HqTU). 

In principle, there is no upper limit to the arrival time 
t arr for the central black hole ejection process. The black 
hole is always there, and on-going star formation (e.g. 
ILu et al.l 12009) provides a constant supply of new stars. 
To derive a typical £ arr , we consider the 'loss cone,' the 
set of orbits which have a distance of closest approach 
within the black hole's tidal radius. For an ensemble of 
stars formed close to the black hole, a few will have orbits 
that interact with the black hole on < 1 Myr timescales 
and so are quickly removed. As a result of dynamical 
interactions with other massive objects or the long-term 
evolution of chaotic orbits within a triaxial potential, the 
remaining stars w ill 'fill' the loss cone with timescales of 
100 Mvr to 1 Gvr (lYu fc Tremaine!l2003t iMerritt fc Poonl 
2004; iWang fc Merrittll2004: iPerets et al.ll2007l) . 

Timescale thus provides an important distinction be- 
tween the central black hole and disk runaway ejection 
processes. The disk runaway scenarios must eject stars 



within the 1-10 Myr lifetimes of massive stars to attain 
unbound velocities. The central black hole can eject un- 
bound stars at any time, however we expect that stars 
formed in the Galactic center will have typical arrival 
times of 0.1-1 Gyr. For the HVSs studied here, an upper 
limit is provided by their finite lifetimes. 

3.3. Comparison with Observations 

From Table [U observed HVSs have i arr » 50-100 Myr. 
The formal error in t arr is likely an underestimate of the 
true error, however. Perhaps the best estimate of sys- 
tematic error comes from comparing the measured stellar 
parameters wit h different sets o f stella r evolution tracks. 
For HVS5, the lEkstrom et "all (|2012t ) tracks for rotat- 
ing stars give a long er ag e of 200 ± 23 Myr, while the 
iGirardi et all (|2002l 120041 1 +0.2 dex super-solar tracks 
give a shorter age of 142 ± 16 Myr. Taking the ±30 Myr 
age spread as the systematic error rules out an ejection in 
the first 10 Myr of HVS5's lifetime at the 3-cr level. This 
confidence level is corroborated by the directly measured 
parameters: the log g and T e g of a 75 Myr old 3.62 M© 
star differ by 2.1-er and 3. 3-cr, respectively, with respect 
to HVS5's present values (see Figure |3|). Thus for HVSS 
we rule out a possible t arr < 10 Myr at the 3-cr level, an 
interesting and important constraint on its origin. 

The hyper-runaways first discovered by iHeber et al.l 
(|2008al) show a contrasting result. HD 271791 is an 
unbound 11 Mq B star. The observed proper motion 
shows it wa s ejected in the dire ction of rotation from the 
outer disk (H eber et al.l l2008ah . The star has an age of 
25 ± 5 Myr and a flight time from the disk of 25 ± 6 Myr. 
Formally, the star has i arr = ± 8 Myr. T he marginally 
unbound 5 Mq B star HIP 60350 is similar (Irrgan g et al.l 
l20ich . It has an age of 45^q Myr and a flight time from 

the disk of 14±3 Myr. Thus, f arr = 31+^ Myr. The short 
arrival times a re consistent with both the supernova ejec- 
tion scenario dPrzybilla et alJl2008aD and the dynamical 
ejection scenario (jGvaramadzd 120091 ) . Contrasting the 
derived arrival times with arrival times for HVSs under- 
scores the usefulness of i arr as a model disc riminant. 

Oth er objects are more ambiguous. iTillich et all 
(2009) discovered the marginally unbound 2.5 Mq A 
star J0136+2425 with a derived age of 245 Myr and a 
flight time of 12 Myr if it comes from the disk. Ac- 
cepting modern Mil ky Way halo mass estimates of ~ 
1.7 x 10 12 M Q (e.g. iGnedin et aLll2010l: iPrzvbilla etaLI 
120101) , it is bound to the Milky Way and may be ex- 
plained as a halo star. The evolved sdB star J1211+1437 
has a flight time that is also a small fraction of its progen- 
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itor's lifetime (which may be many Gyr) (|Tillich et al.l 
2011). Given the ±140 km s _1 uncertainty in the space 
motion, this sdB star is also consistent with being bound 
and thus a normal halo star. 

The HVSs are significantly unbound based on radial 
velocity alone. The four HVSs with known evolutionary 
state discussed here have t arl = 50-100 Myr, times both 
larger than known hyper-runaways and larger than the 
maximum t arr expected in the mechanisms for producing 
hyper-runaways. However, their i arr are close to the ar- 
rival times expected for dynamical interactions with the 
black hole at the Galactic center. 

4. CONCLUSION 

We describe Keck HIRES spectroscopy of HVS5, one 
of the fastest known HVSs with a minimum Galactic 
rest frame velocity of +663 ± 3 km s^ 1 . The obser- 
vations reveal that HVS5 has a projected rotation of 
vsini= 133 ± 7 km s _1 and is thus a main sequence B 
star. Comparing the measured T c g and log g with stellar 
evolution tracks indicates that HVS5 is a 3.62±0.11 M , 
170 ± 17 Myr old star. Given its present distance and ra- 
dial velocity, we calculate that HVS5's arrival time, the 
time between its formation and subsequent ejection, is 
t sn = 105 ± 18(stat)±30(sys) Myr. 

This timescale provides an interesting new constraint 
on the origin of unbound runaways and HVSs. Runaway 
B stars near the disk have t alI — 0-30 Myr, consistent 
with disk ejection scenarios involving a supernova in a bi- 
nary system or a dynamical event among several massive 



stars. The set of B-type HVSs with known evolutionary 
states, on the other hand, have t alI = 50-100 Myr. This 
timescale is difficult to reconcile with any ejection mecha- 
nism requiring a massive star to attain unbound ejection. 
The central black hole ejection scenario, however, allows 
for any i arr . Thus, the derived arrival times for HVSs 
support the black hole ejection model. 

Future progress requires obtaining high resolution ob- 
servations of other HVSs to constrain their age and dis- 
tance. The age distribution of HVSs has important impli- 
cations for the epochs of s tar formation and th e growth 
of the central black hole (jBromlev et al.l 120121 ) . Com- 
bined with future proper motion measurements, we hope 
to directly constrain the full space velocity and place of 
origin of the HVSs. 
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